The complexity of energy eigenstates as a mechanism for equihbration 



Llui's Masanes/ Augusto J. Roncaglia/ and Antonio Acm^'^ 

ICFO-Institut de Ciencies Fotdniques, Mediterranean Technology Park, 08860 Castelldefels (Barcelona), Spain 
^ ICREA-Institucio Catalana de Recerca i Estudis AvanQats, Llms Companys 23, 08010 Barcelona, Spam 

Understanding the meciianisms responsible for tiie equilibration of isolated quantum many-body 
systems is a longstanding open problem. Some of them have been identified, but a complete picture 
is still missing. In this work we obtain a link between equilibration and the complexity of the 
hamiltonian's eigenvectors. We see that hamiltonains with generic eigenvectors equilibrate, and 
provide a simple closed formula for the equilibration time-scale. This allows us to obtain the 
equilibration time-scale for random hamiltonians. We also attempt to quantify the complexity 
of the energy eigenvectors, by considering hamiltonians whose diagonalizing unitary is a quantum 
circuit. We show that when the circuit size is quadratic or larger equilibration is expected, and 
when it is linear or smaller equilibration is not expected. 



INTRODUCTION 

A physical system that has been sitting for a while is 
often described by a thermal or Gibbs state. This pre- 
supposes that, whatever the initial state was, the system 
evolved into a stationary state. But this is impossible for 
closed systems evolving unitarily, unless the initial state 
was already stationary. However, often, the reduced den- 
sity matrix of a subsystem does evolve to a stationary 
state, and stays close to it for most of the time — this 
is called local equilibration. Identifying which conditions 
are responsible for this process is a long-standing ques- 
tion in Physics, both for classical [l|, Q and quantum 
systems Recently, significant advances in the un- 

derstanding of this problem have been achieved due to 
three factors. First, powerful numerical techniques have 
enabled the dynamical simulation of large many-body 
systems @, [3| . Second, the use of quantum- information 
ideas, and in particular of entanglement theory, has pro- 
vided new perspectives into this question [8l4lC| . Finally, 
experiments with ultra-cold atoms have allowed the ma- 
nipulation and observation of many-body systems, with 



high control Interestingly, these experiments have 

challenged the current understanding of these questions, 
since no thermalization was observed in non-integrablc 
systems 13, 13] ■ 

Local equilibration can be explained by dephasing 0, 
under the condition of no degenerate en- 
ergy gaps (no resonant transitions) [l^l- This mech- 
anism could be the complete answer, but it turns out 
that it docs not explain all equilibration processes. Sys- 
tems of quasi-free bosons have infinitely-many degener- 
ate energy gaps, and some of them enjoy local equilibra- 
tion [qI. [itI. . 

In this work, we consider a different mechanism for lo- 
cal equilibration, which is complementary to dephasing 
since it is based on the structure of the energy eigen- 
vectors instead of the eigenvalues. This mechanism is 
independent of the energy spectrum, hence it allows to 
construct hamiltonians which enjoy local equilibration 



despite having many degenerate energy gaps. Addition- 
ally, this mechanism gives reasonable equilibration time 
scales, while the condition of no-degenerate energy gaps 
alone gives unrealistically large values [1, [l^ . 

Our main goal is to explore the relationship between 
the equilibration properties of a system and the com- 
plexity of its energy eigenvectors. First, we show that 
hamiltonians with generic eigenvectors enjoy local equili- 
bration, independently of their spectra and whether they 
have degenerate energy gaps. Second, we show that the 
dynamics of the convergence to equilibrium and its time 
scale do depend on the spectrum, in particular, on the 
Fourier transform of the level density. We calculate this 
function for a random matrix and for the spectrum of an 
integrable system, obtaining their corresponding equili- 
bration time-scales. Third, we attempt to quantify the 
complexity of the eigenvectors of a hamiltonian and re- 
late this to its equilibration properties. In order to do so, 
we consider all hamiltonians whose diagonalizing unitary 
is a quantum circuit with a fixed number of gates. We 
show that, for almost all circuits with number of gates 
quadratic in the system's size the associated hamilton- 
ain enjoys local equilibration. Also, for almost all cir- 
cuits with less than linear number of gates the associated 
hamiltonian does not equilibrate. 



EQUILIBRATION IN HAMILTONIANS WITH 
GENERIC EIGENSTATES 

We consider a system of TV qubits (spin-i systems) 
with hamiltonian 
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where d = 2^. The hamiltonian characterizes the dy- 
namics of the system: if p is the state of the N qubits 
at time t = then p{t) = e~'*^/9e'*^ is the state at 
time t. The diagonalizing unitary U maps a local basis 
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\n) = |ni, . . . , un) to the energy eigenbasis \^n) = U]^), 
for all n S {0, 1}''^. We do not impose any constraint 
on the spectrum of the hamiltonian, and in particular, 
it can have many degenerate energy gaps M. The 
maximal degeneracy of the hamiltonian is denoted by 
g ~ max„ \{n' : En' = En}\. We will show that, in 
hamiltonians with generic eigenstates, the dynamics of 
the convergence towards equilibrium is characterized by 
the Fourier transform of the level density IJ,{E), 
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where the second equality holds for finite-dimensional 
systems. 

Suppose that we are interested in a subset of the N 
qubits — we refer to it as the subsystem, while the rest of 
qubits are referred to as the environment. The Hilbert 
space factorizes as "H = Hs^T-Le, with the corresponding 
dimensions satisfying d = dsds- If p is a state of the N 
qubits, the reduced state of the subsystem is ps = tr_Ep. 
Local equilibration happens when the subsystem evolves 
towards a particular state, and stays close to it for most 
of the time. If the state ps{t) equilibrates to the state 
Ps then, this must be the time-averaged one [lo| 



PS = lim ^ 



T 

dtps it) • 
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As in [sl-fioj , we quantify the departure from equilibrium 
at time t by the trace distance ||ps(t) — Ps||i, which is 
directly related to the probability of distinguishing ps (t) 
from Ps with the optimal measurement |20l . l2l| . The 



following result (proven in Appendix A) bounds the de- 
parture from equilibrium in hamiltonians with generic 
eigenstates. 

Result 1. For any A''-qubit initial state p, almost all 
hamiltonians with a given spectrum {Ei, . . . , Ed} satisfy 
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for all t. 



The meaning of "almost all" is controlled by the free 
parameter e e (0,1), which is an upper bound for the 
fraction of hamiltonians that violate the bound accord- 
ing to the Haar measure [22|. For example, if we set 



e = 0.01 then 99% of the hamiltonians satisfy the above 
bound [i^ . Note that this bound is independent of the 
initial state p, but the set of hamiltonians which violate 
it could depend on p. 

Let us discuss the significance of the three terms in- 
side the brackets of (|4]). The first term depends on the 
spectrum of H and the time t. At time t = the bound 
is useless, since fl{0) = 1. For sufficiently long times, it 
is expected that the phases in the sum ^ cancel each 



other, resulting in a small number, and implying equili- 
bration. Actually, Result 2 below shows that this is the 
case for most of the times, independently of the existence 
of degenerate energy gaps. However, in finite systems 
(d < oo), there are some very special (and very long) 
times ij-cq for which p.{t^cq) ~ m(0)- These are the quasi- 
recurrences, in which the system goes back to a non- 
equilibrium state. In the thermodynamic limit quasi- 
recurrences tend to disappear. In Section „ the quantity 
fl{t) is calculated for some meaningful spectra. The sec- 
ond term also depends on the spectrum of H, and im- 
plies that hamiltonians with high degeneracy cannot be 
warranted to equilibrate. The third term implies that 
equilibration needs the environment to be much larger 
than the subsystem. This condition is necessary in all 
approaches to equilibration known to the authors 0, S~ 

Although quasi-recurrences take the subsystem out of 
equilibrium, the following result (proven in Appendix A) 
shows that, in some circumstances, the subsystem is close 
to the stationary state ps most of the time. 

Result 2. For any A^-qubit initial state p, almost all 
hamiltonians with a given spectrum {Ei, . . . , Ed} satisfy 
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This shows that, in the reasonable regime where g,ds 
ds equilibration is expected for all hamiltonians except 
for a fraction e G (0,1) according to the Haar measure. 
This establishes the existence of a mechanism for equi- 
libration which is only based on the genericness of the 
diagonalizing unitary U, or equivalently, the genericness 
of the energy eigenstates. This mechanism does not rely 
on any condition for the energy spectrum — even in the 
presence of many degenerate energy gaps equilibration 
happens. 



TIME SCALES 

Result 1 shows that the dynamics of the convergence 
to equilibrium is given by the function |/i(<)|, which only 
depends on the spectrum. In this section we calculate 
|/i(i)| for two spectra, each representing an extreme case. 



Spectrum of a random hamiltonian 

Let us consider random hamiltonians sampled from the 
gaussian unitary ensemble According to this, each 



matrix element Hi- 



is an independent random 



variable with probability density 
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It is shown in (2J] that, according to this measure, the 
eigenvalues of H are statisticahy independent from the 
eigenvectors of H, and the diagonahzing-unitary U fol- 
lows the Haar measure [22|, as in results 1 and 2. 

The convergence function corresponding to the spec- 
trum of a random matrix is, in the large-ci limit, 
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where Ji is the Bessel function of first kind (see Appendix 
B). This gives an equilibration time-scale 
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where E-^gy^ is the largest eigenvalue of the hamiltonian. 
Note that both, in terms of N and in terms of E'max, 
the equilibration time ([7]) is much smaller than the one 
obtained for the spectrum of an integrable system, see 
Eq. below. This is expected, since the spectrum 
associated to ([7]) is maximally chaotic. 



Spectrum of an integrable system 

Let us consider hamiltonians with the diagonalizing 
unitary U being generic, and the spectrum being the one 
of the Ising model in a transverse magnetic field h. (Note 
that the hamiltonain of the Ising model need not be in 
this category.) The eigenenergies are parametrized by 
the vectors n = (rti, . . . ,nAr), where e {0, 1} are the 
occupation numbers of the energy eigcnmodes: 
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En = ^nkUj{2TTk/N), 
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k=l 



uj{4>) = y'(/i-cos0)2-hsin2<^ . (9) 
In the i ^ 1 regime we obtain (see Appendix B) 



This gives an equilibration time-scale 
1 1 



(11) 



where -Emax is the largest eigenvalue of the hamiltonian 
(see Appendix B). 



EQUILIBRATION AND COMPLEXITY 

In the theory of quantum computation every algorithm 
can be represented by a quantum circuit, a sequence of 
one and two-qubit gates, which processes the input data 



to generate the output [20(1. This is analogous to clas- 
sical computation, where algorithms can be represented 
by circuits of logical gates. These gates can be seen as 
elementary computational steps, so that harder compu- 
tations require more gates. Following this idea, one can 
quantify the complexity of a unitary matrix by the size 
of the smallest circuit which approximates it to a suffi- 
ciently good accuracy. It seems natural to quantify the 
complexity of the eigenvectors of a hamiltonian by the 
circuit complexity of the diagonalizing unitary U. 



Large complexity 

A central part in the proof of Result 1 consists of per- 
forming an average over all possible unitaries U. How- 
ever, the averaged expression only contains a fourth 
power oi U ® U . Then, if instead of all unitaries one 
averages over a subset being a 4-design, the same result 
is obtained. In [2^ 26 1 strong evidence is provided to 



the fact that random circuits constitute good approxi- 
mations to 4-designs. Assuming this result to be true, 
we can prove prove the following. 

Result 3. Suppose that a given initial state p evolves 
under a hamiltonian H ^ U diag(i?i , . . . , Ed)U\ where U 
is any circuit with C gates. For almost all such circuits 
we have 



\ps{t)~Ps\\ 
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for all t. 



The meaning of "almost all" is controlled by the free pa- 
rameter e S (0, 1), which is an upper bound for the frac- 
tion of circuits with C gates that violate the bound [27 1 . 
The constant a depends on the universal gate set, and it 



is calculated in [25|. Compared to Result 1 there is an 



extra term inside the brackets, which disappears in the 
limit C — > oo. When C is quadratic in the number of 
qubits (or larger). 



C > a'N^ 



(13) 



for a' > 3/a, the extra term is exponentially small in N. 
It is also easy to check that the time average for T — oo 
gives a bound similar to Result 2. Therefore, there seems 
to be a relationship between the equilibration properties 
of the hamiltonian and the complexity of its eigenvectors. 
Note that quadratic scaling is the minimum needed to 
warrant that in the process of diagonalizing all qubits 
"talk" to each other. 
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Small complexity 



CONCLUSIONS 



Consider a hamiltonian with no interaction between 
subsystem and environment: H = Hs®Ie + Is®He 
where Is{Ie) is the identity matrix for the subsys- 
tem(environment). In this case, the reduced density ma- 
trix ps{t) = e^'*^'5p5 e'*^'' does not converge to any- 
thing, unless it is in a stationary state from the begin- 
ning [ps^Hs] — 0. Therefore, interaction is necessary 
for equilibration. The condition of no degenerate energy 
gaps implies that there is interaction across all possible 
bipartitions subsystem-environment It also implies 
local equilibration, independently of the complexity of 
the hamiltonian. Therefore, in order to investigate the 
lack of equilibration, we have to restrict to Hamiltonians 
with many degenerate energy gaps. 

For simplicity, let us consider hamiltonians whose spec- 
trum is that of quasi-free fermions 
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where n = (ni, . . . ,n7v), rik € {0, 1} are the occupation 
numbers of the energy eigenmodes, and Uk are the cor- 
responding excitation energies. Let the subsystem be an 
M-qubit subset of the N qubits, and the environment the 
remaining N ~ M qubits. In the case C = 0, the eigen- 
basis is local |^'„) = |ni, . . . , rijv), each qubit evolves in- 
dependently, the subsystem does not interact with the 
environment, and hence, there is no equilibration. Let 
us see that this is still the case when the complexity is 
sufficiently small. 

Let us compute the probability that a random circuit 
for iV qubits has no gates involving any of M fixed qubits. 
The random circuit is generated by repeating the follow- 
ing process C times: uniformly pick a gate from the uni- 
versal gate set; if this is a single-qubit gate apply it to 
a qubit chosen uniformly from the TV qubits; if this is a 
two-qubit gate apply it to a pair of qubits chosen uni- 
formly. The probability that no gate is applied to any of 
the M qubits satisfies 



In this work we have addressed the problem of equi- 
libration in isolated quantum many-body systems evolv- 
ing under unitary dynamics. We have pointed out the 
existence of a mechanism for local equilibration which is 
based on the complexity of the energy eigenvectors. Un- 
der the action of this mechanism the equilibration time 
scale is given by a nice closed expression, which we could 
explicitly compute for random hamiltonians. 

The relation between equilibration and the complexity 
of solving the dynamics of a physical system emerging 
from our results resembles the situation in classical me- 
chanics, where the notion of integrability plays an impor- 
tant role (H. In fact, there exists a link between equili- 
bration (formalized by weak mixing 0)) and the difficulty 
of solving the dynamics of a system: classical integrable 
systems violate weak mixing, while sufficiently chaotic 
systems satisfy it. Now, if the circuit size is interpreted 
as the complexity of solving the dynamics of a quantum 
system, the resulting picture resembles what happens in 
classical mechanics. Hence, our results may also con- 
tribute to the problem of finding a definition for quan- 
tum integrability — a proposal in terms of computational 
complexity can be found in (28l |. 

Dephasing under the condition of no degenerate energy 
gaps d, and the complexity of the energy eigenvec- 
tors are two independent mechanisms that explain the 
phenomenon of local equilibration. Are there other mech- 
anisms, apart from these two? Does any of them play a 
dominant role in natural phenomena? 

Note added after completion of this work: results re- 
lated to the ones presented here have been obtained in- 
dependently in References [2^ and (soj . 
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In this event, there is no interaction, and hence, no equi- 
libration. Suppose the complexity is sublinear: C < 
with < < 1. If we fix the size of the subsystem Af , in 
the large- limit we have p « 1 — 2M / N'^~'^ , and then 

Result 4. For almost all circuits with sublinear com- 
plexity, the associated hamiltonian with spectrum (|14p 
does not enjoy local equilibration. 



[1] H. B. Callen, Thermodynamics and an introduction to 
thermostatistics (John Wiley & Sons, 1985). 

[2] V.I. Arnold, A. Avez; Ergodic Problems of Classical Me- 
chanics; New York: W.A. Benjamin (1968). 

[3] E. Schrodinger, Annalen der Physik 388, 956 (1927). 

[4] J. Von Neumann, Zeitschrift fur Physik A 57, 3070 
(1929). 

[5] J. Gemmer, M. Michel and G. Mahler, Quantum Ther- 
modynamics (Berlin, Springer, 2004). 

[6] M. Rigol, V. Dunjko and M. Olshanii, Nature 452, 854 
(2008); M. Rigol, Phys. Rev. Lett. 103, 100403 (2009); 



5 



A. C. Cassidy, C. W. Clark, and M. Rigol, Phys. Rev. 
Lett. 106, 140405 (2011). 

M. C. Banuls, J. I. Cirac, and M. B. Hastings, Phys. Rev. 
Lett. 106, 050405 (2011). 

N. Linden, S. Popescu, A. J. Short, A. Winter, Phys. 
Rev. E 79, 061103 (2009). 

M. Cramer, J. Eisert; New J. Phys. 12, 055020 (2010). 
C. Gogohn, M. P. Mueller, J. Eisert; Phys. Rev. Lett. 
106, 040401 (2011). 

L Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

T. Kinoshita, T. Wenger, and D. S. Weiss, Nature (Lon- 
don) 440, 900 (2006). 

S. Hofferberth, I. Lesanovsky, B. Fischer, T. Schumm and 
J. Schmiedmayer, Nature 449, 324-327 (2007). 
S. Goldstein, J. L. Lebowitz, C. Mastrodonato, R. Tu- 
mulka, N. Zanghi, Proc. R. Soc. A 466, 3203-3224 (2010). 
P. Reimann; Phys. Rev. Lett. 101, 190403 (2008). 
A hamiltonian has no degenerate gap whenever its spec- 
trum is such that if Em — En — E^' ^ E„i then m = n 
and m' = n' , or m = m' and n = n' . 
M. Cramer, C. M. Dawson, J. Eisert, T. J. Osborne; 
Phys. Rev. Lett. 100, 030602 (2008). 
A. R. Usha Devi, A. K. Rajagopal; Phys. Rev. E 80, 
011136 (2009). 

Z.-X. Gong, L.-M. Duan; Comment on "Foundation 
of Statistical Mechanics under Experimentally Realistic 
Conditions"; larXiv:1109 .4696 

M. A. Nielsen, L L. Chuang; Quantum Information 
and Quantum Computation, Cambridge University Press, 
Cambridge, 2000. 

The trace-norm of a matrix A is = trV A^A. The 

probability of discriminating between pi and p2 is p = 

\ + illpi - P2||i. 

A. Haar, Ann. Math. 34, 147-169 (1933). 
If the spectrum is degenerate different unitaries can give 
rise to the same hamiltonian. However, the proportion of 
equivalence classes which violate the bound is identical 
to the proportion of unitaries which violate the bound. 
M. L. Metha, Random Matrices, 2nd ed. (Academic 
Press, 1990). 

W. G. Brown, L. Viola; Phys. Rev. Lett. 104, 250501 
(2010). 

F. G. S. L. Brandao, A. W. Harrow, M. Horodecki, to 
appear. 

Recall that the number of circuits with C gates is finite, 
hence there is no need to define a measure over them. 
Also, different circuits can give rise to the same unitary 
or hamiltonian, hence, the bound on the proportion e 
cannot be directly translated to the set of hamiltonains. 
J.-S. Caux and J. Mossel, J. Stat. Mech. (2011) P02023; 
,arXiv:1012.3587. 

Vinayak, M. Znidaric; Subsystem's dynamics under ran- 
dom Hamiltonian evolution, arXiv:1107.6035 
F. G. S. L. Brandao, P. Cwikliriski, M. Horodecki, 
P. Horodecki, J. Korbicz, M. Mozrzymas; Conver- 
gence to equilibrium under a random Hamiltonian, 
larXiv:1108.2985] 



Equilibration bounds 



Proof of Result 1 



Define the linear map 

^M= 5] e"(^"-^"')trB(M'„pM'„0, (16) 

and note that ^t[p\ — Ps{t) ~ PS- The sum 
runs over all pairs of eigenstates n,n' e {1, . . . ,d} with 
different energies E„ ^ En' ■ Any d x d matrix B satisfies 
\\B\\\ < dii{B^B), hence we have the bound 



\ps(t)-Ps\\l<dsiTs{^t[p?)- 



(17) 



Let {|s); s ~ 1, . . . ,^5} be an orthonormal basis of "Hs, 
and {|e); e = 1, . . . , ds} be an orthonormal basis of He- 
Next we consider the case where the initial state is pure 
?/) — \tp){il^\. Some calculation shows 



trs(^^tM') 
■ (se|«'„^/'^'„'|s'e)(s'e'|«'fc?/'*fe'|se') 



(18) 



s.e.s' ,e' 



= {Mt\U (g)U(g)U(g>U(g)U(g)U(g)U(g) t/|right(t)) 
where U is the complex-conjugate of U, and 

(left| = {se,ip, s'e' jipjip, s'e,ip, se'\ 

s 

|right(i)) 



S.e.s' .e' 



En^E„, Ek^E^, 

\n, n' , k, k' , n ,n' ,k , k'). 



22| is denoted by {F{U))u- It is shown in [251 



The average of a f/-dependent matrix F{U) over the Haar 
measure 
that 



where the index tt runs over the 4! permutations 
of four elements, and the normalized vectors |$7r) = 
(K- ^ -f)|</')i5|'/')26|0)37|</')48 are constructed with the 
maximally-entangled state ~ o'^^^^ Sn=i 

between the i"^ and j*^ factor spaces, and the unitary 
which permutes the four left-most factors spaces, and / 
which is the identity in the four right-most factor spaces. 
In the following we will calculate ||/Os(i) — PSII2 exactly, 
but we will bound each of the terms in the expression to 
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illustrate its main features. Some calculus shows 

(left|$,) = d-2(se,V',sV,^|14|V,s'e,V,se') 

s.e,s' ,e' 

($,|right(t)) ^ e"(^"-^"'+^--^'=') 

■ d~^{n,n',k,k'\V^\n,n',k,k'). 

The constraints En 7^ E^' and E^ 7^ E^./ 
imply that almost all permutations tt satisfy 
(n, n', fc, fc'|V^|n, n', fc, fc') = 0. The only permuta- 
tions TT for which (left|$^)(<i>T|right(<)) ^ are analyzed 
in what follows. For the identity tt = (1234) we have 

(left|$(i234)> = d-Hi-s{tTl\i^M)<d~^, 
($(1234)1^8-^^) = E E d-2e'*(^"-^"'+^-^-). 

Ek^E^, 

In a similar fashion 

(lcft|$(3412)) = d-'-tYE{tTl\^M)<d-^ 

($(34i2)|right(t)) = ^ d-2e"(2i=„- 



(lcft|*(432i)) = d-^ {se\se'){s'e'\s'e)^dE\ 



s.e.s .e 



($(432i)|right(t)) = Y d-''<l, 

E„i^E„, 

(lcft|$(i432)) = d-Hl-EitTlli^M) < d-^ 

($(i432)|right(t)) = Y d-^e'*(^"+^-^^"'). 

E„^E„,^Ek 

(lcft|$(32i4)) = d-2^rs(trl;|V')(V^|)<d-2^ 
($(32i4)|rightW) = Y d-2e'*(2^^"+^"'-^'='). 

(lcft|$(i324)) = (rffjd)"' 

($(i324)|right(0) = Y rf^'e"(^"-^'='', 

and for tt = (4231) we obtain the same re- 
sults as for TT = (1324). Note that for 
TT = (3412), (4321), (1432), (3214), (1324), (4231) we have 
|(left|$^)($^|right(t))| < The fact that |||left)|| = d 
and |l|right(t))|| < d^ implies |(left|A|right(0) | < d^. 
Using these two facts, the chain of equalities p8)) . and 
inequality (jl9p . we obtain 



(trs(r!t[V']'))c, < (left|$(i234)|right(i)) 
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ds 



•3\C 



d-'X 



Let p be the not-necessarily-pure initial state (p = 
'^iPi^'i where each tpi is pure). Any real- valued random 



variable X satisfies {X) < {X'^y^'^. Using the above, the 
triangular inequality, and (jl7p . we obtain 



\ps{t) - psW 



< Yp^lmmi 



< Yp^i'^strsm^^r) 



1/2 



< ((left I $(1234) |right(t)) + A + d3^c^ (20) 

Define w = {J2e =e ^^^'^ note that w < gd. Direct 
calculation shows 



(left|$(i234)|right(i)> 

< ^ ^ ^-4 git(B„ 
Ek^E^, 

= \fiit)f + w^d-^ ~2wd-^\tl{t)\'^ 



(21) 



which can be fed in (|20p giving 

|ps(0 - Psili , 



< (lA(t)l 



(7 

4 9^ 6 ,^ „ < 
d2 ^ dE 



d^2- — 



1/2 



We conclude the proof of Result 1 with a simple prob- 
abilistic argument. Let X be a random variable taking 
positive values such that {X) < Xq. If e probjX > a;} 
then (1 — e)0 + ea; < (X), therefore x < x^/e. 



Proof of Result 2 

By applying the previous simple probabilistic argu- 
ment to (|20|) , we obtain that almost all hamiltonians with 
complexity C satisfy 

\\ps{t)--ps\\i (22) 
< ^(^(ieft|ci>(i234)|right(0) +^^2-^ + A' 

Note that 



lim ^ j dt(left|$(i234)|right(t)) 



(23) 



Y Y d'^SiEn - En' + Ek -Ek') < I , 



where S{En — En' + E^ — Ek') is a Kronecker delta. This 
bound, inequality (j22p . and the convexity of the square 



7 



root, imply 



lim - / dt\\ps{t) - psWi 

T^oo T Jq 
wl/2 



< 



lim - / rfi(lcft|$(i234)|riglit(i)) 



4 \i/2 



^■^2'— + — 



e V dE dE 



1/2 



(24) 



This shows Result 2. 



Proof of Result 3 

Result 3 ean be proven by following the same steps 
as in the proof of Result 1, but replacing the average 
over the Haar measure by the aveage over cireuits with 
complexity C. Let /C(C) be the set of A^-qubit cireuits 
with C gates, from a particular universal gate set. The 
average of a function F{U) over U G IC{C) is denoted by 



{F{U)) 



ueK.(c) 



It is shown in [251 that 



u®u®u®u®u®u®u®u) 



< ^|<i>,)($,| + A^A, 



(25) 



where A = 1 — a In 2 /N and the constant a > depends 
on the universal gate set. The hermitian matrix A satis- 
fies ll^lloo < 1 and A\^Tr) = for all tt. The range of tt 
and the vectors |$7r) are defined in the proof of Result 1. 



erated by the probability distribution P(H_) correspond- 
ing to the gaussian unitary ensemble [2^. According 
to this, each matrix element Hij is an independent ran- 
dom variable; the elements in the diagonal Ha £ R have 
probability density P{Hii) = 7r~^/^e~-^"; the elements 
not in the diagonal Hij = Hji G C have probability den- 
sity P{Hij) = (2/7r) e~^l^'jl^. This is equivalent to say 
that the diagonalizing-unitary U oi H follows the uniform 
distribution over unitaries (the Haar measure (22j). and 
independently, the spectrum of H follows the probability 
density 



PiE,,...,Ea) = ae-^'=^^^ Y[ \E, ~ Ej\^ (26) 

l<i<j"<d 

where Ei G (— oo, oo), and a is a normalization constant. 
That is, eigenvalues and eigenvectors are independent 
random variables. 

Since, what appears in expression is \p{t)\ to the 
fourth power, we are going to calculate the average 



(ImWI^)h = jdEi---dEdPiEu...,Ed)\mf- 

A standard trick within random matrix theory is that, 
with probability almost one, the value of |/i(i)|"' for a 
randomly chosen spectrum is very close to the above av- 
erage. 

For the next, it is useful to define the n-point 
correlation- function 

dl f 

i?„(i?l, En) = y- TT / rfi?„+i • ■ ■ dEd P{Ei, . . . , Ed) 

[d-n)\ J 

where < n < d |24|. The sum 



Calculation of |/i(t)| 
Spectrum of a random matrix 

In this section we calculate the convergence function 
|/i(<)| for the spectrum of a random hamiltonian i?, gen- 



ijkl 



H 



can be split in the four terms where i,j,k,l are: (i) all 
different, (ii) two of them equal; (iii) three of them equal 
or two pairs equal, and (iv) all equal. These four terms 
arc separated in 



d''{\m\')H 



JdEi dE2dE3dE4 R4{Ei,E2,E3,E4) c'^i^^-E^+E^-E^) 

dEi dE2dE3 R3iEuE2,E3) ^e''^^E^-E,-E,) ^ ^-iti2E^~E,-Es) ^ ^^itiE,-E,) 
dE,dE2 R2{Ei,E2) \c^t2(E^-E,) ^ ^^.tiE^-E,)-]^ ^ [2d^ - d] . 

I 



(27) 



It is shown in |24| that the n-point correlation functions can be written as 

Rn{Ei, ...,E„)= dct[K{E„E,)],^,=i_^^ (28) 
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where 



(29) 



fe=0 



(|3T1) is bounded by d. Since (1^(01^)// is equal to ((27)) 
divided by d^, all terms are of order l/d or smaller, ex- 
cept for the single term with a four- fold product of ([50]) . 
This implies that in the large-d limit we have 



where (pk{x) are the eigenfunctions of the quantum har- 
monic oscillator. To see how the determinant works in 
consider the example 

R2iEuE2) - K{Ei,Ei) K{E2,E2) - K{Ei,E2)^ . 



When substituting 
objects of the form 



{\m\')H 



2Ji{tV2d) 



\ W2d J 
Spectrum of the Ising model 



(36) 



in (P7)) . we obtain products of 
dE K{E, E) e'*-^ (30) 



The eigenenergies are parametrized by the vectors 
n = (rii, . . . ,npf), where Uk € {0, 1} are the occupation 
numbers of the energy eigenmodes: 



N 



and 

JdE, ■ ■ ■ dE,, K{E,,E2) ■ ■ ■ K{E,,,E,) eHti^^i+-+*"S„ 



En = ^nuuj{2^k/N), 

k=l 

'^(4') — y cos 0)^ -|- sin^ 



(37) 
(38) 



= tr[P 



•Pe 



(31) 



for n = 2, 3, 4 — where P = X)fc=o \ is the projector 

onto the d-dimensional lower-energy subspace of the har- 
monic oscillator, and X is the position operator. For any 
pair of bounded operators A, B, the Cauchy- Schwartz in- 
equality tells 



Then 



-2N 



N 

n 

fc=i 

N 



k/N) 



2"^ (1 + cos[tuj{2T:k/N)]) 



\tiiA^B] \ < y^tr[AtA] tr[BtB] . (32) 
This can be used to bound ([51]) for the case n = A: 
tr [(P e'*i^P e'*^^) (P e"-^^P e'*^^)] 



fc=i 



2-^exp 



2ir 



— /d(/)ln(l + cos[iw((?!))]) 
27r 



(39) 



where the approximation holds for large N . Let us ana- 
lyze this expression in the small-i and large-t limits. For 
t ^ 1 we have, up to fourth-order terms. 



< ^tr [e-i*2Xp e-"i^P2 e^i^P e"^-^] tr [P e-"3^p e'*^^] 

(33) 



4' 



< 



tr [P] tr [P] J tr [P] tr [P] < d. 



By setting t4 = and ^3 ~ 0, we obtain the same bound 
for n = 3, 2. It is shown in 2j] that in the large-d limit 
we have 



\a{l + cos[tu{<j))]) wln2- (l + h- 2h cos\(t)\) — 
which gives 

|/i(t)|2«e-*'^(i+"')/^ 
For t 1 wc have 



(40) 



A^ln 1 



1 



This is called "the semi-circle law". In this limit, the 
integral ([5(1)) can be evaluated 



27r JO 

2-^exp[iVln(l)] = 2"^ 



2ir 

d(j) cos[ta;((/)) 



(41) 



fdEK{E^E)c^'^ « d'Hl^^ < d, 
J tV2d 



(35) 



where the inequality follows from the convexity of the log- 
arithm, and the approximation holds when the integrand 
oscillates heavily, that is when t ^ 1. In the large- iV 
limit, the largest eigenvalue is 



where Ji is the Bessel function of the first kind. When 
substituting (|28)) in (|27| . we obtain several terms, each 
being a product of objects of the form (|30)) or (|3T)) . Ac- 
cording to (1551) and (|35p . each of these factors (1501) or 



Erria^c 



N 
2^ 



[h — cos ( 



sm 



(42) 



which is proportional to N . 



